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I.  INTRODUCTION 


Th«  purpose  of  this  report  is  to  deteraine  the  validity  of  Rockwell 
Collins'  assumptions  that  g-sensltlve  terms  MU^2  and  MDb2  sre  small,  stable,' 
and  negligible.  This  report  contains  the  results  of  tests  and  evaluations 
conducted  by  the  Research,  Development,  and  Engineering  Center  to  confirm  the 
magnitudes  of  these  terms.  Appendix  A  contains  a  detailed  discussion  of  tests 
and  test  results. 

II.  BACKGROUND 

Rockwell  Collins  assumes  that  the  magnitude  of  the  gyro  acceleration  sen¬ 
sitive  terms  MU^  and  MUB2  in  equations  (1)  through  (6)  below  are  small, 
stable,  and  negligible.  It  is  well  known  that  the  GSA2  and  GSB2  acceleration 
sensitive  terms  in  these  equations  are  highly  temperature  dependent. 

The  Collins  Phase  III  Multi function  Sensor  Inertial  Measurement  Unit 
(IMU)  is  used  in  the  evaluation.  Prelaunch  self-calibration  measurements  are 
made  at  three  stationary  orientations.  These  orientations  for  gyro  2,  at  each 
measurement  position,  are  shown  in  Figure  1.  Figure  1  also  provides  the  com¬ 
ponents  of  acceleration  (G^,  G2,  G3).  The  measurement  equations  for  each  gyro 
output  data  axis  are  as  follows: 

The  measurement  equations  for  Collect  1  orientation  are: 


RA2^  ■  B^2  “  GSA2 
RB  2  *  B|2  ^  GSA2 


Rotation  of  the  multisensor  to 
measurement  equations: 

RA2^  *  Bj^j  +  GSA2 

RB2^  •  B|2  +  GSA2 


Rotation  of  the  multi sensor  to 
measurement  equations: 


3 

RA2  ■  8^2  +  GSA2 
3 

RB2  ■  B32  +  GSA2 


*  Gt  +  GSB2  *  G2  +  MUb2  *  G3  -  (1) 

*  G2  +  GSB2  *  Gt  -  MU^2  *  G3  +  W2  (2) 

Collect  2  position  yields  the  following  static 

*  Gi  +  GSB2  *  G2  -  MUb2  *  G3  +  Wx  (3) 

*  G2  -  GSB2  *  Gx  ♦  MU^2  *  G3  +  W2  (4) 

Collect  3  position  yields  the  following  static 

*  G3  +  GSB2  *  G2  +  MUb2  *  Gi  +  W3  (5) 

*  G2  -  GSB2  *  G3  -  MU^2  *  Gt  +  W2  (6) 


The  six  measurement  equations  above  are  used  to  solve  for  g— sensl tlve 
(6SA2,  GSB2,  M0^2 •  MDB2)  and  bias  (B^2,  BB3)  terms*  The  on-axis  g-sensltlve 
(GSA3)  term  is  obtained  by  subtracting  equation  (1)  from  equation  (3). 

GSA2  *  Gx  -  (RA22  -  RA21)  +  2  +  MDb2  *  G3  -  Wx  (7) 

The  cross-axis  g-sensltlve  (GSB2)  term  is  obtained  by  subtracting 
equation  (4)  from  equation  (2): 

GSB2  *  Gj  -  (RB21  -  RB22)  +  2  +  MD^  *  G3  (8) 

Observe  that  the  stability  of  the  GSA2  and  GSB2  terms  depend  on  the  va¬ 
lidity  of  the  MUB2  and  MtJ^  assumptions  and  how  well  the  acceleration  com¬ 
ponents  are  defined.  For  this  evaluation,  the  magnitudes  of  the  G2  and  G3 
gravity  components  are  reduced  to  approximately  zero. 

Equations  (9)  and  (10),  derived  from  equations  (1)  through  (4),  are  used 
to  compute  the  gyro  bias  (B): 

B^  *  (RA22  +  RAj1)  *  2  -  GSB2  *  G2  (9) 

bb2  -  (RB22  +  RB2l)  »2-  GSA2  *  G2  -  W2  (10) 


Note  the  bias  stability  Is  dependent  on  the  predictability  of  the  GSA2,  GSB2 , 
and  G2  terms. 

To  calculate  sensor  2  MU^  and  M0B2  terms,  sensor  2  gyro  equations  were 
set  up  in  the  following  matrix  form: 
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The  following  augmented  matrix  was  obtained,  row  reduction  procedures 
applied,  and  the  solution  to  the  six  gyro  parameters  found. 

The  sensor  2  augmented  matrix  is 
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where 

cY  -  ra2  (Cl)  +  wx 

c2  -  RB2  (Cl)  ~  w2 
C3  -  RA2  (C2)  -  Wi 
C4  -  W>2  (C2)  -  W2 
C5  -  RA2  (C3)  -  W3 
C6  -  RB2  (C3)  -  W2 


The  sensor  2  matrix  gyro  parameter  solutions  for  MU&2  and  MUjj2  are  aa  follows 


and 


MOb2  -  (2SA23  -  RA2L  -  RA22  -  2W3)  Gl  +  (RA2L  -  RA22  +  2WX)  G3  (12) 

2(0^  +  G32) 

Ill.  EVALUATION  OF  THE  MU^  AND  M0b2  ACCELERATION  TERMS 

During  chese  evaluations,  an  environmental  chamber,  placed  over  the 
multifunction  IMU,  provided  the  mechanism  for  controlling  the  ambient  tem¬ 
perature  environment.  Twenty-four  prelaunch  self-calibration  tests  were  made 
with  the  instrument  mounted  In  the  navigation  orientation.  Three  self¬ 
calibration  tests  were  made  at  each  ambient  temperature  setting.  Forty-five 
minutes  between  each  self-calibration  tests  were  allowed  for  the  instrument  to 
thermally  restabilize.  The  sequence  was  repeated  at  the  eight  designated  tem¬ 
perature  settings. 

The  test  results,  showing  the  computed  mean  magnitude  values  and  one 
sigma  randomness  of  the  gyro  bias,  MU^2,  and  MUb2  drifts,  are  summarized  In 
Table  1.  Figures  2  and  3  show  line  graphs  of  the  bias  data  plotted  for  the 
eight  temperature  settings.  The  matrix  line  graphs  (M)  (Figs.  2  and  3)  repre¬ 
sent  the  bias  drift  obtained  by  using  equations  (9)  and  (10).  In  utilizing 
equations  (9)  and  (10),  the  MU a2  and  MUB2  acceleration-sensitive  terms  cancel 
and  the  bias  Is  dependent  only  on  the  predictability  of  the  GSA2  *  G2  and 
GSB2  *  G2  terms.  Line  graphs,  derived  from  equations  (5)  and  (6)  which  assume 
the  magnitudes  of  MU^  and  MUB2  terms  are  zero,  are  plotted  also  In  Figures  2 
and  3  with  the  appropriate  matrix  line  graph. 

The  line  graphs  (Figs.  4  and  5)  depict  the  magnitudes  of  the  MU^2  and 
MUb2  acceleration  sensitive  terms  derived  by  using  equations  (11)  and  (12). 
These  line  graphs  represent  the  terms  chat  Collins  assumed  were  small,  stable 
and  negligible. 
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Figure  2 .  Bias  A  drift  va  temperature 


Figure  3.  Bias  B  drift  vs  teaperature 


MUa  drift  vs  teaperature 


IV.  CONCLUSIONS 


Fro*  Che  tests  and  evaluations  conducted,  data  were  obtained  (Table  1) 
which  confirms  that  the  MU^  and  MUB2  acceleration  sensitive  terns  are  tem¬ 
perature  dependent.  The  magnitude  of  the  MU^  term  varies  from  23.8  to 
-37.2* /hr/g  with  a  maximum  uncertainty  of  10.3* /hr/g  one-sigma.  The  magnitude 
of  the  MDg  term  varies  fro*  3.6  to  -24.3* /hr/g  with  a  maximum  uncertainty  of 
6.4* /hr/g  one-sigma.  In  conclusion,  the  bias  drift  data  also  reflects  the 
Impact  of  assuming  that  MD^  *ud  MUB2  terms  are  equal  to  zero.  The  line 
graphs  (Figs.  2  and  3)  clearly  illustrate  significant  differences  in  the  bias 
drift  magnitudes  when  the  terms  MU^  and  MUB2  are  assuned  small,  stable,  and 
negligible  when  they  are  not. 

V.  RECOMMENDATIONS 

The  MU^2  and  MDB2  acceleration/temperature  sensitive  drift  rates  do  play 
a  significant  role  in  computing  the  gyro  bias  if  equations  (5)  and  (6)  are 
used.  To  eliminate  these  terms,  it  is  recommended  that  bias  equations  (9)  and 
(10)  be  used. 

In  using  equations  (9)  and  (10),  the  computed  bias  values  are  still  depen¬ 
dent  on  the  predictability  of  the  following  terms:  GSA2  and  GSB2  accelera¬ 
tion/temperature  sensitive  drifts,  the  acceleration  components  defined  by  the 
multi sensor  coordinate  frame,  and  the  components  of  earth  rate  for  the  defined 
coordinate  frame.  The  ability  to  accurately  predict  these  values  will  deter¬ 
mine  the  performance  limitations  of  the  multi sensors.  It  is  recommended  that 
the  values  of  these  terms  be  well  established  and  fitted  to  polynomials  to 
determine  if  the  multlsensor  performance  can  be  significantly  Improved  by 
software  modeling  these  dependent  variables. 
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GLOSSARY 


Symbol 

®A2»  ®B2 
Gl»  ®2»  c3 

GSA2  *  Gt 

GSB2  *  Gt 
*®A2  *  G1 

*®B2  *  G1 

ra2,  rb2 

Wi  (i  -  1,  2,  3) 

i  (i  -  1,  2,  3) 


Definitions 

Gyro  drift  bias  for  axes  A2  and  B2. 

Components  of  acceleration  corresponding  to  the  coordinate 
frame  defined  by  the  multi sensor  2. 

Gyro  g-sensitive  drift  due  to  acceleration  along  the  angu¬ 
lar  rate  axis. 

Cross-axis  g-sensitive  drift. 

Gyro  g-sensltlve  drift  due  to  multlsensor  A-axis  accelera¬ 
tion  sensitivity. 

Gyro  g-sensltlve  drift  due  to  multlsensor  B-axls  accelera¬ 
tion  sensitivity. 

Plckoff  output  for  gyro  axes  A2,  B2. 

Components  of  earth  rate  corresponding  to  coordinate  frame 
defined  for  multi sensor  2. 

Superscripts  1,  2,  3  represent  Collect  positions  1,  2,  3. 
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APPENDIX 

DETAILED  DISCUSSION 


Th«  Collins  Phase  III  multifunction  sensor  uses  a  self-calibration  scheme 
to  calculate  gyro/accelerometer  parameters.  To  accomplish  this,  the  two  sen¬ 
sors  which  make  up  the  multifunction  sensor  are  rotated  to  three  orientations. 
Figure  1  shows  sensor  2  rotation  positions.  At  each  position,  raw  drift  rate 
data  Is  collected  for  axes  A2  and  B2.  The  drift  rate  data  Is  composed  of  g- 
sensitlve  terms,  bias  terms,  and  earth  rotation  terms.  It  la  necessary  to 
determine  the  values  of  these  gyro  parameters.  For  a  particular  latitude,  the 
earth  rotation  terms  are  known.  Therefore,  the  unknowns  reduce  to  the  g- 
sensltlve  parameters  and  the  bias  parameters. 

For  each  sensor,  at  each  orientation  position,  two  gyro  measurement 
equations  can  be  written  In  terms  of  the  g-sensitive  parameters,  the  bias 
parameters,  and  the  earth  rate  components.  Sensor  2  gyro  measurement  equa¬ 
tions  for  the  three  collect  positions  are  written  as: 

(COLLECT  POSITION  1) 

RA21  -  MUb2  *  g3  +  *A2  ”  GSA2  *  Gl  +  GSB2  *  G2  -  Wx 
RB21  -  -MD^  *  G3  +  Bb2  +  GSA2  *  G2  +  GSB2  *  Gj^  +  W2 
(COLLECT  POSITION  2) 


RA22  -  -MUB2  *  G3  +  ♦  GSA2  *  Gi  +  GSB2  *  G2  +  Wx  (A-3) 

RB22  -  MUA2  *  G3  +  832  +  GSA2  *  G2  -  GSB2  *  G^  +  W2  (A-4) 

(COLLECT  POSITION  3) 

RA23  -  MUb2  *  +  GSA2  *  G3  +  GSB2  *  G2  +  W3  (A-5) 

RB23  -  “MGa2  *  G1  +  B^  +  GSA2  *  G2  -  GSB2  *  G3  +  W2  (A-6) 


At  this  point,  there  are  six  equations  and  six  unknowns.  The  unknowns  are 
HUb2>  MUa2  »  GSA2,  GSB2,  Ba2  and  Bb2.  Collins  makes  the  assumption  that 
MUb2  and  MU^  are  small  and  therefore  negligible.  This  assumption  reduces  the 
number  of  unknowns  from  six  to  four  with  the  number  of  equations  staying  at 
six.  Therefore,  sensor  2  gyro  equations  for  the  three  collect  positions 
reduce  to: 


(A-l) 

(A-2) 


(COLLECT  POSITION  1) 


RA21  -  Ba2  *  6SA2  *  Gi  +  GSB2  *  G2  -  Wi  (A-7) 

M21  -  Bb2  +  GSA2  *  G2  +  GSB2  *  Gi  +  W2  (A-8) 

(COLLECT  POSITION  2) 

BA22  -  Ba2  +  GSA2  *  Gi  +  GSB2  *  G2  +  Wi  (A-9) 

RB22  -  Bb2  +  GSA2  *  G2  -  GSB2  *  Gi  +  W2  (A-10) 

(COLLECT  POSITION  3) 

BA23  -  Ba2  +  GSA2  *  G3  +  GSB2  *  G2  +  W3  (A-ll) 

RB23  -  BB2  +  GSA2  *  G2  -  GSB2  *  G3  +  W2  (A-12) 


Making  the  assumption  that  MUg2  and  MDx2  are  zero  reduces  the  unknowns  to 
GSA2,  GSB2,  Ba2  and  Bb2*  At  this  point,  there  are  six  measurement  equations 
and  four  unknowns. 

The  following  procedure  was  used  by  Collins  to  solve  the  four  unknown 
gyro  parameters  for  sensor  2.  First,  Collins  subtracts  equation  (A-7)  from 
equation  (A-9)  and  solves  for  GSA2 •  Next,  Collins  subtracts  equation  (A-10) 
from  equation  (A-8)  and  solves  for  GSB2*  Collins  then  substitutes  the  pre¬ 
viously  solved  parameters  for  GSA2  and  GSB2  into  equations  (A-ll)  and  (A-12) 
and  solves  for  Ba?  and  Bb2»  respectively.  Collins'  solutions  for  these  param¬ 
eters  are  written  as  follows: 


GSA2  -  RA22  -  RA21  -  2Wi 

(A-13) 

2G1 

GSB2  -  RB21  -  RB22 

(A-14) 

2G1 

Ba2  "  RA23  -  GSA2  *  G3  -  GSB2  *  G2  -  W3 

(A-15) 

Bfi2  "  RB23  -  GSA2  *  G2  +  GSB2  *  G3  -  W2 

(A-16) 

A-2 

$>;cr:iv-/^%xv:vs»s^v:v^ 


If  MDa2  and  MDb2  are  not  zero,  Collins'  solution  equations  (A-13)  through 
(A-16)  would  change  as  follows: 


GSA2 

-  ra22  - 

•  RA21  -  2Wi  +  MDB2 

*  G3 

(A-17) 

2G1 

GSB2  -  RB21  -  RB22  +  M0a2  *  G3 

(A-18) 

2G1  G1 

®A2 

-  RA23  - 

GSA2  * 

G3  -  GSB2  *  G2  -  W3 

-  MDg2 

*  Gi 

(A-19) 

®B2 

-  RB23  - 

GSA2  * 

G2  +  GSB2  *  G3  -  W2 

+  M&A2 

*  Gi 

(A-20) 

Solution  equations  (A-17)  through  (A-20)  present  a  problem.  Gyro  parameters 
GSA2  and  GSB2  are  now  a  function  of  known  parameters  and  unknown  parameters 
(M0a2»  MDb2)«  If  Collins'  assumption  of  M0b2  and  MDx2  equals  zero  is  Incor¬ 
rect,  then  the  solution  scheme  Is  invalid.  In  laboratory  testing,  G2  and 
G3  are  very  small  while  Gi  is  in  the  +lg  field.  Therefore,  if  MUg2  and  MD^2 
are  not  zero,  the  only  terms  that  would  appear  incorrect  are  the  bias  terms 
Ba2  and  Bg2. 

At  this  point,  it  was  decided  to  check  the  assumption  that  MGb2  and  MUa2 
equal  zero.  If  the  assumption  was  proven  Invalid,  it  was  decided  to  show  the 
effect  on  the  gyro  parameters. 

Before  the  assumption  is  made,  there  are  six  equations  and  six  unknowns. 
Mathematically,  It  is  possible  to  calculate  the  unknown  parameters.  Therefore, 
sensor  2  gyro  equations  were  set  up  in  matrix  form  (p.  3).  Obtaining  the 
augmented  matrix  (p.  4),  and  applying  row  reduction  procedures,  the  solutions 
to  the  six  gyro  parameters  were  found.  The  matrix  gyro  parameter  solutions 
for  sensor  2  are: 


M0a2  “  (RB21  +  RB22  “  2RB23)  Gi  +  (RB22  -  RB2*)  G3  (A-21) 

2(G^2  +  G3 2) 


MUb2  -  (2RA23  -  RA21  -  RA22  -  2W3)  Gi  +  (RA^  -  RA22  +  2^)  G3  (A-22) 

2(Gi2  +  G32) 


GSA2 


RA22  “  RA21  -  2Wi  +  MUB2  *  G3 
2G,  Gi 


(A-23) 


A-3 


.  *.  * 


(A-24) 


GSB2  -  RB21  -  RB22  +  MDa2  *  G3 
2Gi  Gl 


Ba2  -  (RA21  +  BA22)  -  GSB2  *  G2  (A-25) 

2 

Bb2  "  (BB21  +  RB22  -  2W2)  -  GSA2  *  G2  (A-26) 

2 


Note  that  each  parameter  is  a  function  of  that  previously  solved  for  parame¬ 
ters  and/or  known  raw  data. 

To  determine  whether  MD\2  and  MUjj2  values  were  small,  self-calibration 
runs  were  conducted  over  a  temperature  range  from  -21.5  °F  to  140  °F.  The  data 
in  Table  1  shows  that  H0^2  varies  from  -37.2  */hr/g  to  23.8  */hr/g  and  MUg2 
varies  from  -24.3  */hr/g  to  3.6  s/hr/g  over  the  temperature  range. 

Table  1  also  shows  the  effects  of  the  Invalid  assumption  on  and  Bg2* 

In  laboratory  testing,  Gi  is  approximately  equal  to  +1G.  To  solve  for  B^2  and 
Bg2»  Collins  uses  equations  (A-15)  and  (A-16),  respectively.  To  include  the 
effects  of  MU&2  and  M0jj2  m  Che  bias  terms,  equations  (A-15)  and  (A-16)  were 
modified  to  form  equations  (A-19)  and  (A-20),  respectively.  Using  values  from 
Table  1,  it  can  be  shown  that  the  modified  solutions  for  Ba%  and  Bg2  agree 
with  the  matrix  solution. 

Collins  uses  Collect  Position  3  equations  (A-ll)  and  (A-12)  to  solve  for 
B^2  and  Bg2«  Collect  Position  2  equations  (A-9)  and  (A-10)  and  Collect 
Position  1  equations  (A-7)  and  (A-8)  could  have  been  used  to  solve  B^2  and 
Bb2*  It  was  decided  to  solve  for  B^2  and  Bg2»  using  equations  (A-7)  through 
(A-10),  and  to  compare  the  results  to  the  results  Collins  calculated  using 
equations  (A-ll)  and  (A-12).  It  was  also  decided  to  modify  equations  (A-7) 
through  (A-12)  to  take  into  account  for  MUa2  and  MUg2«  1°  simplify  the  com¬ 
parison,  only  the  results  for  B&2  will  be  shown.  B^2  was  solved  for  using  the 
following  equations: 

Proa  matrix  solution  equation  (A-25) 

BA2  -  (RA21  +  RA22)  “  GSB2  *  G2  .  (A-27) 

2 


Prom  equation  (A-7),  Collect  1 

Ba 2  -  RA21  +  GSA2  *  Gi  -  GSB2  *  G2  +  Wi.  (A-28) 


i 


From  equation  (A-9),  Collect  2 


BA2  “  *A22  “  GSA2  *  Gi  -  GSB2  *  G2  -  Wi.  (A-29) 

Proa  Equation  (A-ll),  (Collins'  solution) 

BA2  "  B^23  “  GSA2  *  G3  -  GSB2  *  G2  -  W3.  (A-30) 

Froa  equation  (A-l ) 

BA2  -  RA21  +  GSA2  *  Gl  -  GSB2  *  G2  +  Wi  -  MDb2  *  G3.  (A-31) 


Proa  equation  (A-3) 


BA2  •  BA22  -  GSA2  *  Gi  -  GSB2  *  G2  -  Wl  +  HDB2  *  G3. 


(A-32) 


Proa  equation  (A-5) 


BA2 


-  GSA2  *  G3  -  GSB2  *  G2  -  W3  -  MDB2  *  Gi. 


(A-33) 


Equation  (A-27)  Is  the  matrix  solution  for  B^2*  Equations  (A-28)  through 
(A-30)  are  solutions  for  B&2  from  Collins'  Collect  positions  1,  2,  3,  equa¬ 
tions.  Equations  (A-31)  through  (A-33)  are  solutions  for  B&2  from  Collins' 
Collect  positions  1,  2,  3,  modified  to  take  Into  account  for  MU^2  and  MDg2* 

In  laboratory  tests,  Gi  Is  approximately  equal  to  1G  while  G2  and  G3  are 
approximately  zero.  Therefore,  equations  (A-31)  and  (A-32)  are  approximately 
the  saae  as  equation  (A-28)  and  (A-29),  respectively. 

Self-calibration  tests  weie  performed  over  a  small  temperature  range. 

The  raw  data  and  gyro  parameters  were  substituted  Into  equations  (A-27) 
through  (A-33).  Table  A-l  shows  the  results  of  these  self-calibration  tuns. 
Notice  that  all  the  results  for  B&2  are  approximately  the  same  except  for 
Ba2  solved  from  equation  (A-30).  This  is  the  equation  Collins  uses  In  their 
solution  scheme.  The  result  for  Ba2  solved  from  equation  (A-33)  is  Collins' 


solution  taking  into  account  MO^  and  MDB2*  This  modified  solution  agrees 
with  the  other  solutions  for  BA? .  The  same  comparisons  can  be  shown  for  Bg2« 

From  the  results  shown  in  this  appendix,  it  is  concluded  that  MU^  and 
MUB2  are  parameters  that  need  to  be  calibrated*  Neglecting  these  parameters 
cause  some  gyro  parameters  to  be  calculated  Incorrectly.  It  is  recommended 
that  the  matrix  solutions  be  used  for  the  gyro  parameters  of  both  sensor  1  and 
sensor  2  of  the  multifunction  sensor. 


BA2  Eq.(A-27)  BA2  Eq.(A-28)  Ba2  Eq.(A-29)  BA2  Eq.(A-30)  Bo  Eq.(A-31)  Bo  Eq.(A-32)  Bo  Eq.(A-33) 
(Matrix  (from  (from  (used  by  (Modified  (Modified  (Modified 


-V-'^V-VJV  •->  -<-g  • 


19.4812  19.4754  19.4865  14.7815  19.3996  19.5623  19.4796 
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